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Abstract: We study the thermodynamics and the small fluctuations of a linear dilaton 
black hole, in an S-dual version of the near-horizon limit of type IIB NS5 branes. The 
thermodynamical analysis shows that the black hole is in a Hagedorn phase, with marginal 
stability. The dynamical analysis confirms that the speed of sound of the dual theory 
vanishes and that there are no evident instabilities at the level of supergravity. We clarify 
the physical meaning of some singularities of the retarded correlator that were thought to 
lead to instabilities. 
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1. Introduction 

Little string theories (LSTs) are non-local theories without gravity that describe the dy- 
namics of NS5 branes in the limit where the string coupling vanishes [1-3]. Although the 
excitations on the brane decouple from the bulk, the theory on the NS5 branes remains 
strongly coupled. LST admits a holographic description in terms of a near-horizon lin- 
ear dilaton geometry [4]. At finite temperature, the only geometries known analytically 
correspond to a Hagedorn phase, so the temperature is fixed at the scale of the LST. 

There are several features that make these theories interesting. The low energy de- 
scription of LST in type IIA is given by the (2,0) superconformal theory, whereas for LST 
in type IIB the low energy description is in terms of Super Yang Mills theory with (1,1) 
supersymmetry. The holographic description of LST generalizes the more usual AdS /CFT 
correspondence [5]. However, the theory has typical properties of a string theory like T- 
duality, and the holographic dual allows to compute off-shell quantities. This is in contrast 
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with critical string theories, where so far it is possible only to calculate the S-matrix. An- 
other interesting point is that the linear dilaton black holes are closer relatives to black 
holes in asymptotically flat space [6-8]. Then, understanding Hagedorn LSTs may be a 
step towards the implementation of holographic ideas to less exotic black holes. 

A key issue is the thermodynamical stability of the system. In a holographic construc- 
tion there is an identification of degrees of freedom on both sides of the duality, so the 
existence of thermodynamical instabilities should be reflected in the dual geometry along 
the lines of the Gubser-Mitra conjecture [9]. Indeed, the conjecture has been confirmed to 
hold for backgrounds with finite negative specific heat, cy < 0, where there is a dynamically 
unstable mode that is associated to an imaginary speed of sound in the dual theory [10] 

2 _ dP _ S_ 

The Hagedorn phase of LSTs is marginally stable, but at weak coupling it can be shown that 
quantum corrections render it unstable [11-14]. At strong coupling, black hole solutions 
that would correspond to LST on compact spaces have been shown to be thermodynami- 
cally unstable with negative specific heat [10,15,16]. In the flat space limit, the specific heat 
diverges, cy — > — oo, so that the speed of sound vanishes and the background is marginally 
stable [10]. Then, according to the Gubser-Mitra conjecture, the geometry should not show 
any instability at the classical level, although it could be unstable if the effect of quantum 
corrections were taken into account. 

The absence of the tachyon associated to the sound mode was confirmed in [17]. How- 
ever, the computation of the holographic Green's function [17-19] reveals an analytic struc- 
ture that seems incompatible with a retarded correlator and maybe introduces other in- 
stabilities. If this were the case, the conjecture would fail. Our purpose in this paper is 
to clarify this issue using a particular setup. We analyze the retarded Green's function 
in the linear dilaton black hole and compute the spectral function above the mass gap for 
the operators dual to the dilaton and the metric. We will show that there are no evident 
instabilities, and explain the origin of the analytic structure found previously. We also find 
a new set of modes describing energy loss that was not discussed previously. 

The outline of the paper is as follows. In section |2| we will introduce a family of 
black holes that are the holographic duals to a LST at different volumes and study their 
thermodynamics. We will concentrate on the infinite volume case and, in section ||, we will 
give the equations of motion for small metric and dilaton fluctuations and describe some 
common qualitative features of their effective potentials. In sections |I| and ||, we will discuss 
the decay of fluctuations in the black hole geometry through quasinormal modes and argue 
that there are no dynamical instabilities. Then, we will relate the previous results to the 
scattering matrix and the holographic computation of two-point functions in section ||. We 
will present our conclusions in section ^. Technical details of the calculations can be found 
in the appendix. 

2. Thermodynamics of the Hagedorn phase 

LSTs are non-gravitational string theories that describe the dynamics of NS5 branes. They 
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have a density of states that grows exponentially with the energy 

n(E) = E a e pHE . 

The Hagedorn temperature Tjj = 1/Ph is a limiting temperature of the system. The 
partition function 

Z = j dEfl(E) e~ pE 

becomes formally divergent beyond this point and a first order phase transition will take 
place, although the final state is not known for LSTs. The first order phase transition 
is a characteristic of the canonical ensemble. In the microcanonical ensemble the energy 
of the system can be increased until the Hagedorn temperature is reached. Then, adding 
further energy does not increase the temperature and it is invested in exciting internal 
degrees of freedom instead. To a first approximation, the temperature does not change 
at all, but, depending on the sign of a, there can be a slight slope that gives positive 
specific heat if a > or negative specific heat if a < 0. Thus, the Hagedorn phase could 
be thermodynamically stable or unstable. It is a highly non-trivial question whether the 
Hagedorn phase of a given theory would be stable or not. 

In order to study the LST at strong coupling we will use a holographic description in 
terms of an S-dual geometry sourced by a configuration of D5 branes. The solution [20,21] 
has been used as the gravitational dual of the theory living on N c five-branes wrapping an 
S 2 . At low energies the theory flows to SU(N C ) N = 1 SYM in flat space [21]. We 
will consider a deformation of this theory, with Nf = 2N C five-branes that introduce 
fundamental degrees of freedom in the dual theory [22] . These geometries were constructed 
in [22] and a linear dilaton black hole solution at the Hagedorn temperature of the dual 
LST was found there as well. As we will discuss later, the presence of the linear dilaton 
implies that the low energy degrees of freedom of the field theory do not really decouple 
from the higher LST modes . Nevertheless, the geometry close to the horizon should still 
be associated to their holographic description. 

In [25], a new class of linear dilaton black holes with compact horizon was found with 
temperature T higher than the Hagedorn temperature in flat space Tu = t. — A , • The 
metric in the Einstein frame is 

ds = ea 



f(r)dt* + tfdnl + + N c a'dYi 

= 2(r -r) 



/(r) = l- e ^o-r J) < p = ( f >0 + r . (2.1) 

where r is the holographic radial coordinate, (j) is the dilaton and I5 is the internal space, 
whose explicit form is unimportant for our analysis. Although quantum corrections will 
start to be important as r is increased, we are interested in the classical behavior and we 
will ignore this. In order to remain in a weakly coupled regime, one should perform an 
S-duality for some value of r and continue along a geometry with decreasing dilaton. This 



is not very important for our analysis, as will become clearer in section 3.1 



*At zero temperature there are special cases where there is total decoupling, as was studied in [23,24]. 
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The horizon is at r = and has a three-sphere S s component of radius R = Nh/2 
proportional to the flux of Ramond-Ramond 3-form on the sphere 

F {3) = N h . (2.2) 

The temperature of these black holes is independent of the radial position of the horizon 



T(fl) = OTV 1 + § = r "V 1 + l >T - (2 - 3) 

We can think of the T = Th black hole solution with flat R 3 horizon as the formal limit 
R — ► oo of the solution (2.1) 



lim T(R) = T+ (2.4) 

H—*oo 

This limit is consistent with the vanishing of the Ramond-Ramond 3-form flux density on 
the three-sphere, 

N h 2 . . 

The entropy of these solutions was computed in [25]. We compute the energy ( |A.6[ ) 
and find that 

E = TS = Ce 2r °V(S 3 )^, C = N f 2<h W V ^) . (2 . 6) 

Th 87rGio 

Here V(S S ) is the volume of the S 3 and V(Y§) the volume of the internal five-dimensional 
space. Our results agree with the computation at infinite volume [26] when the R — > oo 
limit is taken. 

We can interpret each black hole solution as a state of fixed entropy and volume in the 
microcanonical ensemble. The temperature is clearly 

If we fix the volume, we can change the energy without changing the temperature by tuning 
the position of the horizon, so 



dE /v 

We can use the last result to find the density of states, 

i = ^|§®^)=SV^. (2.9) 
Then, all these backgrounds are marginally stable, a = 0. 

3. Fluctuations in the black hole 

Small perturbations around the black hole can be studied using the linearized equations 
of motion in the background. This allows to study the hydrodynamics of the dual field 
theory using the methods pioneered in [27] and, more generally, to evaluate the spectrum 
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of quasinormal modes in the black hole background, which describe dissipation in the dual 
theory [28,29]. Hydrodynamic modes in dual LST backgrounds were studied in [17]. 
The metric of the flavored Nf = 2N C planar black hole is the R —* oo limit of ( |2.l| ) 



ds 2 = e5 



-f{r)dt 2 + dx 2 + N ca'j^ + N c a'dY£ 



(3.1) 



where t, x are the coordinates of Minkowsky space. 

For simplicity, we will work in units where a' = 1. We also make the change of variables 
r = — (logu)/2, so the black hole factor simplifies to f(u) = 1 — u. Notice that r — > oo 
corresponds to u — > 0. 

The metric fluctuations split into three independent channels classified by the SO(2) 2 
unbroken group of rotations in space: scalar, vector and tensor. Scalar perturbations mix 
with the dilaton. 

The equations of motion for plane-waves in four-dimensional space are derived in Ap- 
pendix [B]. Let us use dimensionless frequency and momentum in units of the temperature 
to = uj/(2itTh), q = qj(2i\Tii)- From the metric and the dilaton equations, we can build 
two diffeomorphism-invariant quantities Z sounc [ and Z^, obeying 

7 „ , q 2 ^ - 2 "> 2 7 , 

sound (l--u)(q 2 (u-2) + 2tr 2 ) sound 



-4u 2 (l - u)q 2 + (u 2 -3u + 2)q 4 + (3u - 4)q 2 ro 2 + 2ro 4 
4u 2 (l-u) 2 (q 2 (u-2) + 2tr> 2 ) 

,4 o„2^,2 



z 



sound 



4(q 4 - 2q 2 »t, 2 ) _ 
3(l-n)(q 2 (n-2) + 2ro 2 ) ^ _U ' ^ > 



7" 1 7> i *> 2 - (1 - u)q 2 ^ 

Z * - T^u Z * + 4u 2 (l-u) 2 Z + ~ ° (3 - 3) 



and 



If one sets to zero, equation ( |3.2| ) becomes the same as the sound mode equation found 
in [17] for six-dimensional LST. When the spatial momentum is zero, q = 0, the two 
equations decouple and the sound equation reduces to the scalar equation 

The same equation is also obtained for the tensor component of the metric. 
For the shear mode, we find the equation 

y„ /A, ^ 7/ { s q 2 (l ~U) ~ /N_ n f**\ 

J.hearW + (1 _ u) (q2 (1 _ u) _ tf) Z * hear W " 4(1 _ u) 2 M 2 ^hearW - . (3.5) 

All these equations coincide with the equations found in [17] for six-dimensional LST. 



"For non-zero momentum. 
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3.1 Schrodinger potentials for classical modes 

In order to understand better the asymptotic conditions on the solutions to the equations 
of motion, it is convenient to write them in a Schrodinger form using a Regge- Wheeler 
tortoise coordinate. The resulting "effective potentials" can also be used to study the 1- 
loop corrections to the black hole free energy due to supergravity modes. This was done 
in [13] for the system of flat D5-branes at weak coupling. In that case, 1-loop corrections 
lift the Hagedorn degeneracy of the black hole. 

Let us take Z^{u) = yjuij}(u) and change variables to the Regge- Wheeler tortoise 
coordinate v = log(u) — log(l — u). In this new coordinate, spatial infinity corresponds to 
v — > — oo, while the black hole horizon is at v — ► oo. Equation (3.2) has now the form of a 



Schrodinger equation (— <9 2 + V(v)*j ip = Eip, where the energy is 



and the potential is 



(3.6) 

v{v) = i-" 2 M + ^(i— M) (3 . 7) 

The potential is positive and monotonic, and it is equal to (1 + q 2 )/4 at infinity and 
vanishes at the horizon. Solutions can be expanded in ingoing and outgoing plane waves at 
the horizon. At infinity there are two possible situations. For E < (l + q 2 )/4 or in 2 — q 2 < 1, 
one solution blows up and the other is exponentially decreasing. On the other hand, for 
E > (1 + q 2 )/4 or It) 2 — q 2 > 1, the two solutions correspond to outgoing or ingoing plane 
waves. 

We can repeat the same analysis for the sound and the shear mode, obtaining similar 
results 



Zsound = y/u(q 2 (u - 2) + 2tU 2 )V> , Zshear = y/u(?0 2 ~ q 2 (l - uj) l\) 

l- u 2 ( v ) + q 2 (l- u ( v )) 2q 2 (2tt, 2 - q 2 )(l - u(v))u 2 (v) 

V soun d{v) = h 



4 {q 2 (u(v) -2) + 2tt) 2 ) 2 

l-u 2 (v) + q 2 (l-u(v)) q 2 (l-u(v))u 2 (v)(q 2 (l-u(v)) + 2tv 2 ) 

VvhearlV) = ; 1 ^ 

K ' 4 4(g 2 (l-^))-tr 2 ) 2 

The potentials have similar asymptotics. The behavior close to the horizon is typical 
of black holes, while the asymptotics at infinity are due to the non-trivial dilaton profile. 
The height of the potential is given by the slope of the linear dilaton background and it is 
independent of the sign, so that the same potential will be found for the S-dual geometry 
as we had anticipated. 

We have computed the equations and the potentials explicitly for the Nf = 2N C 
background, but other black hole backgrounds asymptotically conformally equivalent to 
flat space and with a linear dilaton will have the same asymptotic behavior of the effective 
potential, flat and vanishing at the horizon and constant at large values of the radial 
coordinate, e.g. [17-19]. 

Therefore, we expect that the following analysis will produce qualitatively similar 
results in more general cases. 
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4. Quasinormal modes and speed of sound 



Localized fluctuations in a black hole background tend to spread and disappear as time 
passes. The energy that initially is concentrated in some finite region falls into the black 
hole or escapes to infinity, so that the local energy density vanishes. This process is 
described by the quasinormal modes (QNMs). In order to define QNMs, suitable boundary 
conditions have to be imposed at the horizon and at spatial infinity. The condition at the 
horizon is that the mode should be infalling. In terms of the effective potentials computed 



in section 3.1, this means that we pick the plane wave moving towards the horizon v — > oo 

-imt+iwv/2 



~ e 



At spatial infinity v — * — oo, we have two possible situations, depending on the four- 
dimensional mass shell tn 2 — q 2 . If to 2 — q 2 < 1, then the modes are below the potential 
barrier, so we must pick the vanishing solution at infinity (Dirichlet boundary condition). 
However, if to 2 — q 2 > 1, then the modes are above the barrier and the right boundary 
condition is to pick the outgoing plane wave. The asymptotic form of the mode is then 

ft 2 - q 2 < 1 ~ e -itTrf+Vl+q 2 -raV2 
tr 2 — q 2 > 1 ~ e -int-W tt> 2 -q 2 -lf/2 

For Imro < 0, the QNMs describe fluctuations that decrease exponentially in time, hence 
their association with dissipation. On the other hand, modes with Imto > would be in 
principle interpreted as instabilities of the background. The quasinormal spectrum of the 
dilaton and shear modes can be obtained analytically, because their equation of motion 
can be reduced to a hyper geometric equation. Although this is not the case for the sound 
mode, we are able to show that the speed of sound vanishes, v s = (4.20). This is in 



agreement with our expectations from the thermodynamical analysis: a vanishing speed of 
sound signals a marginal (in)stability. 

4.1 Solutions of the hypergeometric equation 

In this section we will establish the notation that will be used in the rest of the paper. 
Consider the hypergeometric equation 

u(l - u)y" + (c-(a + b+ l)u)y' - aby = , (4.1) 

The local solutions close to the singular points are given in terms of hypergeometric func- 
tions iF\ (a, b, c; x) = 1 + ^x + 0{x 2 ), x^O. At spatial infinity (it 0+), 

y(u) = a ^ 0) (n) + b yft(u) , 
Vo \u) = 2 Fi{a,b,c;u) , 

2/5 1} (n) = u 1 - 2 F!{a + 1 - c,b + 1 - c,2 - c;u) . (4.2) 

At the horizon (u — * 1~), the local solutions are 

y{u) = aiyf\u) + bxyf\u) . 
(u) = 2 -Pi (a, b, 1 — c + a + b; 1 — u) , 
yW(«) = (l-u) c - a - b 2 F l (c-a,c-b,l + c-a-b;l-u). (4.3) 
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Both sets of solutions are related by the connection coefficients 



(4.4) 



whose explicit analytic form is 

Coo 
Coi 
Cio 



r(l-c + a + b)r(l-c) 
T(l -c + 6)r(l -c + a) 
T(l -c + a + b)T(c - 1) 

rpj 

r(l + c-q-6)r(l-c) 

r(i - a )r(i - b) 



r(i + c -a-b)r( c -i) 
r(c — ojr(c — a) 

When the values of o, b and c are such that a connection coefficient vanishes, then one 
of the hyper geometric solutions becomes a polynomial. On the other hand, for the values 
where the connection coefficient has a pole, the hypergeometric series diverges. 

4.2 Dilaton and tensor mode 

The equation of the dilaton with nonzero momentum, (|3.3| ), reduces to a hypergeometric 
equation via 

Z4u)=u c / 2 (l- U )- im / 2 y(u) (4.6) 
where the parameters of the hypergeometric equation are 



c = 1 + y/l + q 2 - tt> 2 (4.7) 

and 

c .tu 1 + y 7 ! + q 2 - tu 2 - ito 
a = b = i — = 4.8 

2 2 2 y 1 



With this choice for the dilaton perturbation fl4.6|) , the infalling boundary condition cor- 
responds to y(u) ~ 1 close to the horizon, u — > 1~, v — ► +oo, that is the yj "* solution 

^-(1 -«)-*" / 2 yS 0) (n)~e^/ 2 

The Dirichlet condition at spatial infinity, u — > + , v — » — oo corresponds to y(u) ~ 1, 

where 

is the normalizable solution, for values Im It) = 0, to 2 — q 2 < 1, corresponding to the principal 
branch of the square root Re yl + q 2 — ro 2 > 0. One can choose to set the branch cut for 
\J\ + q 2 - m 2 at ImlD = 0, ro 2 > q 2 + f, with Re \f\ + q 2 - ro 2 > on the sheet that 
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first sheet 



+ 

^^^^ 



second sheet 



+ 



+ 



+ 



Figure 1: Branch cuts of yl + q 2 — in 2 on the complex tr plane. Branch points are at 
to = ± \J\ + q 2 , the first sheet corresponds to Re y/l + q 2 — tt> 2 > and the second sheet to 
Re v 7 ! + q 2 - w 2 < 0. The signs above and below the branch cuts correspond to the sign of 
Im \Jl + q 2 - tr 2 . 



includes the physical frequency domain and Re yl + q 2 — tr 2 < on the second sheet. 
This choice implies that next to the branch cuts, 



sign(Im yl + q 2 — ro 2 ) = =F sign(Re to), Im to — > d 



(4.9) 



on the first sheet and the opposite signs on the second sheet. The sign of the imaginary 
part of the square root for each sheet is represented in figure |l[ 

The solution is initially denned on a real line above the branch cuts Imrrj —* + . From 
( |4.9| ) this implies that V< also corresponds to the outgoing solution, for ro 2 > q 2 + 1. 
Indeed, when Re to > y 7 ! + q 2 , 



-ivot 



~ e 



-iRerot+i Imy'l+q 2 -fD 2 v/2 



(4.10) 
oo and it is 

(0) 



so the solution corresponds to a plane wave moving towards infinity v 
straightforward to see that the same holds for Rero < 0. 

Therefore, we find that quasinormal modes correspond to solutions where y^' is pro- 
portional to j/q . The condition is that the connection coefficient Cqi = 0. Using (4.5) and 

(El) 



i + Vi + q 2 -tt) 2 



-n, n 



0,1,2, 



(4.11) 



There are no solutions to this equation on the physical sheet. This agrees with previous 
results found in related Little String backgrounds [17]. 

It is possible to do an analytic extension of the solution to the second sheet. This is 
equivalent to a different choice of solution at spatial infinity 



~ u 



In this case quasinormal modes correspond to solutions where yj "* is proportional to ■ 
The condition is that the connection coefficient Cnn vanishes 



1 — c + a 



1 - \J\ + q 2 - ro 2 - ito 



-n, n = 0, 1, 2, . . . 



(4.12) 
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This condition is denned now on the first sheet of the square root, notice that it can be 
obtained from ( 4.11| ) by flipping the sign of the square root. In this case there is an infinite 



discrete set of solutions 

" = - 2 ' WTT) + W+UJ ' " = °' h % ■ ■ ■ ■ (4 ' 13) 

For q = the n = mode is actually not in the spectrum (Coo 7^ 0) and the full set of 
quasinormal modes have Imtu < 0. When q ^ 0, the n = mode appears above the real 
axis to ~ +iq 2 /2 and as we increase the momentum, more modes move to the upper half 
plane. 

4.3 Shear mode 

We have given an equation for the shear mode ([T^) in terms of the gauge-invariant quantity 
Z shear ■ F° r the analysis of the quasinormal modes we will work instead with non- invariant 
quantities H' tx and H' zx . The location of the quasinormal modes is itself gauge-invariant, 
so it is irrelevant which variable we are using. As a consistency check, we will see that H' tx 
and H' zx have the same quasinormal spectrum. The equations for the shear variables are 

tvH' tx (u) + q(l - u)H' zx {u) = , 

jt« ( \ H L(u) c\xoH tx {u)+to 2 H zx {u) 

zx \ ) (i_ u ) + 4(l-u) 2 u 2 

ttii ( \ q 2 H tx (u) + qtoH zx {u) 
4(1 — u)u A 

We can solve for H' zx from the first equation and for H zx from the third. Deriving the third 
equation with respect to u and plugging back the expressions for H zx and H' zx , we find 
a hypergeometric equation for H' tx . Similarly, we can solve for H' tx in the first equation 
and for H tx in the second. Then, deriving the second equation with respect to u and 
substituting the expressions of H tx and H' tx , we find the same hypergeometric equation for 
H' zx . This shows that the quasinormal spectrum is the same for both. 
The parameters of the hypergeometric equation are 

c= 1 + Vi + q 2 -tt 2 , (4-15) 

and 



c .to -1 + y/l + q 2 - tn 2 - iro 

2~ ~ Z 2"~ 2 



2 2 2 v ; 

The analysis follows in the same way as for the dilaton. The quasinormal modes are 
determined by the conditions 

a = —n, n = 0, 1, 2, . . . 

b = -n, n = 0,1,2, ... (4.17) 
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The solution of a = is a mode with ro = — iq 2 /2, which is precisely the dispersion relation 
of a hydrodynamic mode in the holographic dual. Apart from this mode, the above system 
of equations yields the same quasinormal spectrum that we found for the dilaton ( 4,13; ), 
Therefore, modes in the upper half plane appear for finite values of the momentum, but 
these modes with Imro > are all on the second sheet. 

4.4 The speed of sound 

Let us set the dilaton to zero. The sound equation is a more complicated equation with 
four regular singularities at u = 0, 1, 2 — 2ro 2 /q 2 , oo. The change Z = u p l 2 (\ — u) a / 2 y(u), 
with p = 1 + \J\ + q 2 — ro 2 and a = — itv gives the Heun equation 

II ( 1 b 6 \ I ~ Q _ n (A Q\ 

V + \u + + u - 2 + 2ro 2 /q 2 J V + u(u - l)(u - 2 + 2ro 2 /q 2 ) y ( } 



where 



7 = p = 1 + y 7 ! + q 2 - K> 2 
5 = cr + l = l — ito 

e = -2 

a = j3 = + ^ + = ^(-1 - ito + \J\ + q 2 - ro 2 ) 

2 

Q = \ ~ f (3» + 2i(l + Vl + q 2 -tf 2 )) 
tr 2 



+ « 2 - v / i + q 2 -ft) 2 + itt)(i + Vi + q 2 -^ 2 )) 



q' 2 



(4.19) 



There is an explicit constant solution y(u) = 1 corresponding to the sound mode. Such 
solution appears when a = f3 = Q = 0, that in this case is satisfied for the dispersion 
relation ^ 

tt> = -zi- . (4.20) 
Therefore, the speed of sound is zero, v s = 0, in the Nf = 2N C holographic dual theory. 



5. Retarded two-point function in the black hole 

We have found a set of QNMs with Imro > for q > 0. In order to clarify the issue of the 
QNMs with positive imaginary part, we will follow a similar analysis to [30]. We consider 
the Fourier transform of solutions of the Schrodinger equations introduced in section 3 for 
times t > 

poo 

ip(W,q,v)= e ivot ii){t,^v)dt, lmro>0. (5.1) 

As in [30] the solution ^(tt^q,^) can be thought of as the evolution by an appropriate 
Green's function of a perturbation that starts at t = 0. The Green's function G(ro, q; v, v') is 
determined by the choice of two linearly independent solutions of the Schrodinger equation, 
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■0< and tp > , that satisfy the appropriate boundary conditions at v = — oo and v = +00 
respectively. The explicit form is 

q . Wj ^ = ^M^<h^«w>) (52) 

Vv (it), q) 

where v < = min(«, v'), v > = max(f , v') and W^tu, q) is the Wronskian of the two solutions, 

d d 

iW(xo,q) = V'<(tt),q,w)^'i/ ; >(n',q,v) - ^^<(*f , q, *#>(ttJ, q, v), (5.3) 

As we saw in section [O], the potential vanishes at v = +00 and asymptotes to (l+q 2 )/4 
at v = —00. Physical solutions tjj (to, (\) should be at least delta-normalizable in the bulk. 
This imposes the following boundary condition as v —* +00: 

Notice that this is the prescription to compute a retarded correlator, since the black hole 
does not radiate classically. Conversely, as v — > —00, we find that solutions can be a linear 
combination of 

For values ImtD = + , ro 2 < q 2 + 1, corresponding to Re sj\ + q 2 — ro 2 > 0, only ifi< is 
normalizable. For values of the frequency ro 2 > q 2 + 1, V>< also corresponds to an outgoing 
plane wave, as was explained in section |||. 

The retarded bulk correlator ( |5.2p for the dilaton and the scalar modes has the following 
form 

G(u, u>) = i (^<">) (c ' 1)/2 [( 1 -^<)( 1 -^>)]"' ro/2 2i?i(a; ft> c . u< ) 2jFl(a , 6j l_ c+a+ 6; l- u> ) , 

(1 - CjGoi 

(5.4) 

where we used the fact that V>< = u^^il-ny^^y^ and ^ = nM^l-uT^/ 2 ^. 

Since V< and V> are solution of a Schrodinger equation, their Wronskian (|5.3f ) is 
constant. This constant can be evaluated by considering the limit u — > 0. Using 

and 

yj 0) = Cooyf + CoiyJ 13 = Coo 2F 1 (a, b, c; u) + On u 1 ^ 2 Fi(a + 1 - c, 5 + 1 - c, 2 - c; u) 
we find 
iVT = u c {l - u) l - in C m 



~ = (1 - c)Qn = - r(1 " C + Q + 6)r(C 



r(a)r(6) 

(5.5) 



which leads to (5.4). 
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The Wronskian vanishes when the connection coefficient between the two solutions, 
V>< and has a zero. The inverse of the Wronskian (for the scalar mode) is proportional 
to 

/ / \ 2 

p / l+Vl+q 2 -w 2 -t nj 



r( v / f + q 2 - w 2 )r(i - ixo) 

An important subtlety in the Green's function are the singularities of the solutions. The 
solutions are proportional to hypergeometric functions whose series coefficients diverge for 
some complex values of ft), implying that there is no solution. This happens when c < or 
l + a + 6 — c<0in the first sheet or 2 — c < or 1 - a — 6 + c < in the second sheet. 
The bad frequencies are 

to = —in, to 2 = 1 — n 2 + q 2 , n = 1, 2, . . . 

This is taken into account by the singularities of the Wronskian, that cancel against any 
such spurious singularities. Notice that those solutions correspond to the poles in the 
overall factor, as can be seen in (15. 1 



W ~T(1- c + a + b)T(c) (5.7) 

The singular solutions are for instance at the origin of singularities on the real axis that 
were found in the holographic computation of the dual theory Green's functions [17-19]. 

Then, the only genuine singularities come from the zeroes of the Wronskian, and one 
can show that with the above choice of solutions and ifi>, the corresponding Green's 
function has no singularities when it is analytically continued on the Im to > half plane 
(4.11). Therefore, this Green's function can be identified with the retarded correlator. 



The full analytic extension of the Green's function includes the Im to < half plane 
and the second Riemann sheet that can be reached through the square root branch cuts 
Imtr = 0, ttJ 2 > q 2 + 1. This amounts to flipping the sign of the square root in the 
Wronskian fl5.6| ), which is equivalent to changing the ip< solution to ip<. 

Singularities on the Im to < half plane on either sheet correspond to solutions of the 
Schrodinger equation that diverge as v — * oo and that can be associated to absorption by 
the black hole at late times, as for the usual QNMs studied in AdS spaces. In this case, all 
the poles are on the second sheet 

n(n + 1) q 2 

to n = -2i-A r+i^^- 1 tt> n = l,2, .... 

(2n + l) (2(2n + l))' ' ' 

Singularities on the upper half plane of the second sheet Imru > 0, correspond to modes 
that diverge as v — > — oo, so they also describe the loss of energy in the bulk but as it 
escapes towards spatial infinity. This is clearly different from QNMs in AdS spaces, where 
the geometry acts effectively as a box and the energy loss is produced only by absorption 
by the black hole. 

In summary, it is possible to define a retarded Green's function, analytic in the upper 
half frequency plane of the physical sheet. The QNMs with ImttJ > are all in the second 
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sheet, so they are related to processes with no black hole absorption and where all the 
energy escapes to infinity. 

For the shear mode a similar analysis leads to 

p / -l+^l+g 2 -m 2 -iw \ p / Z+^l+^-xo 2 -! 



W{to, q)- 1 ^ / V '- . (5.8) 

lVl + q 2 -W 2 )r(l-ittO 

The spectrum of QNMs is the same on the second sheet. On the first sheet, there is a 
single pole associated to the shear diffusion in the dual theory 

■ q 2 

6. Scattered waves and holographic spectral function 

We now understand the physical meaning of the unusual (as compared to AdS) QNMs 
from the perspective of the linear dilaton geometry. We now want to interpret them in 
terms of the holographic theory. The simplest approach is to use the relation between the 
greybody factors for the fields in the black hole geometry and the spectral function of the 
dual theory [31,32]. 

When the black hole is at equilibrium with a thermal gas, the emitted radiation is the 
same as the absorbed radiation, so the greybody factor is simply the transmission coef- 
ficient. Then, we have to solve the scattering problem of radiation coming from spatial 
infinity. The scattering problem in the context of zero temperature linear dilaton back- 
grounds was studied previously in [6], where a relation between the scattering matrix and 
holographic Green's functions was proposed. In this case, there are real pole singularities 
that are associated to normalizable states and that should be extracted from the scat- 
tering amplitude in order to compute the Green's function, following a LSZ reduction as 
the one proposed in [33]. Normalizable states were localized at the bottom of the linear 
dilaton geometry, that in our case has been swallowed by the black hole. These states have 
disappeared from the spectrum and the reduction is not necessary in our case. 

For to > \/ q 2 + 1, radiation enters the black hole from spatial infinity and is partially 
reflected and partially transmitted 

v -> oo , e ivov l 2 

where Af n = Coo an d ^4^ ut = Coi are the connection coefficients of the hypergeometric 
equation given in the section ||]. 

For negative frequencies ro < — \/ q 2 + 1, the scattered waves should describe the time 
reversal of the absorption process 

v -> oo , e - inv / 2 

.--CO, A- e iV^^lv/2 +A7 -i^^iv/2 ( 6 - 2 ) 
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where A- = C\\ and A 



out 



c 



10- 



Comparing the 'probability current' J = —i(i/j*d v i/j — c.c. ) at the boundary and at the 
horizon, we can define the transmission and reflection coefficients. Prom the expression of 
the connection coefficients (f4.5|), and using the relations 



r(i + »x)r(i 



TTX 



IX) 



sinh(7ra;) ' 



rii 



IX 



Til 



IX 



7T 



COsh(7Tx) 



r(l + x) = xT{x) 



it turns out that the coefficients for scalar and shear perturbations are the same 



± 



T 



± 



K 



±tt> 



M ± I 2 
\At\ 2 



cosh 2 


_f0 2 -q 2 


-1=F») 


cosh 2 


_|0 2 -q 2 


-lift)) 



l\A 



±12 



1 -TV 



± sinh(7rro) sinh(7TY / tt' 2 — q 2 — 1) 



cosh 



l±to] 



(6.3) 



(6.4) 



Perturbations close to the mass-shell of the Hagedorn temperature tr 2 — q 2 ~ 1 suffer a 
very small absorption by the black hole ~ 1, ~ 0. Higher mass perturbations on 
the other hand are easily absorbed, since IZ^ 1 — > and — > 1, with a large mass limit 
(to — > ±oo) 



1 



cosh (7Ttt)) 



0, T ± ~ tanh 2 (7rro) -> 1 . 



(6.5) 



Using the expression of the transmission coefficients T =t (tt)), we can write down the 
spectral function for the dual theory. Assuming time reversal invariance, the spectral 
function should be an odd function of the frequency, p(vo) = —p(—tv), and it should be 
positive for positive frequency. Then, 

p LST (tf ) oc r + (w)G(tt) - Vi + q 2 ) - r-(to)e(-tt) - Vi + q 2 ) , (6.6) 

where 0(x) is the unit step function. Notice that the spectral function vanishes below the 
mass gap and that there is no discontinuity since r ± (± v / TTq 2 ) = 0. This follows from 
our analysis, and otherwise the spectral function would be ill-defined, since the factors 
are oscillatory below the mass gap. Notice also that the zeroes coincide with the values 
where there are no solutions in the geometry. Finally, although the explicit form of ( |6.4| ) 
has the information about the values of the QNMs that describe energy loss towards spatial 
infinity, their interpretation as resonances is unclear, since the spectral function does not 
admit a simple analytic extension outside the real line. 

If the gap is sent to zero, the spectral function becomes smooth 

PLSriyo) ~ tanh 2 (tu) , 

while for a very large gap M or a large to, it is almost like a step function 

Plst{k) ~ ®(w — M) — 6(-tt) - M) , 

recovering the zero temperature behavior. So, in the Hagedorn phase there is a general 
suppression, especially of the states close to the zero momentum mass gap. 
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The QNMs that contribute to this spectral function correspond to those in the second 
sheet of the retarded Green's function in the black hole. Only the hydrodynamic modes 
are not captured by the little string theory degrees of freedom. They can be seen in 
computations of the retarded Green's function of the field theory using the boundary 
action for modes vanishing at spatial infinity, as in [17-19]. However, the Green's function 
computed in this way presents spurious singularities on the real and imaginary axis. This 
suggests that the usual holographic prescription for the computation of correlators should 
be modified in this case. 

7. Discussion 

The linear dilaton black hole shows a quite peculiar analytic structure in its correlators due 
to the mass gap. In other holographic setups like AdS/CFT, the retarded Green's functions 
only have singularities in the lower half frequency plane. Exceptions to this rule when a 
quasinormal frequency crosses the real axis to the upper half of the complex plane translate 
into dynamical instabilities. This is the case of the instabilities appearing for near-critical 
embeddings of flavor D7 branes in AdS§ as the quark mass is varied [29,34,35], the unstable 
sound mode of backgrounds above the Hagedorn temperature [10] or instabilities of the 
Gregory-Laflamme kind [36]. 

The reason is that while a quasinormal mode usually has a divergent behavior at spatial 
infinity, when it crosses the real axis it becomes a bounded state in the geometry. It has 
finite energy and thus belongs to the physical spectrum, but the energy is negative, which 
means that it corresponds to an instability. In the linear dilaton geometry, the modes that 
cross the real axis always blow up at spatial infinity, so they never belong to the physical 
spectrum and do not produce any instability. A technical way to see this is that they are 
on a 'second sheet' of the retarded Green's function. Other apparent instabilities, including 
poles on the real axis, are actually not present. They correspond to special values of the 
frequency where there are no well-defined solutions to the equations of motion. 

This solves the apparent inconsistency of the Gubser-Mitra conjecture for holographic 
constructions with the results found in [17-19]. The marginal stability of LST at the 
Hagedorn temperature is confirmed by the dynamical stability of the dual linear dilaton 
geometry. Of course, quantum corrections can drastically change this result, but this is 
beyond the realm of classical supergravity, which was our main interest. The holographic 
computation of the spectral function above the mass gap confirms this picture, and no 
instability appears. 

The discrepancy between the bulk analysis and the holographic computation of Green's 
function using the boundary action suggests that the latter needs to be modified. A possible 
way would be to use the bulk Green's function to construct a properly defined bulk-to- 
boundary propagator. The holographic Green's function could then be found from the 
convolution of two such propagators in the bulk or by finding the boundary-to-boundary 
limit. In the presence of the black hole, it may be necessary to perform an analytic extension 
beyond the horizon, as in [38]. 
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A. Energy of T > T H black holes 

The energy of the black hole solutions can be evaluated using the results of [37]. Consider 
a constant time hypersurface St in the black hole background. The metric on St is given 
by 

dr 2 



cfe Is, = ea 



R 2 dU 2 + 7— + N c a'dY 2 



N c a'R 2 . . 



Take Sf° to be a surface of constant and large radial coordinate value r = r max in St, 

'■bh 



which we can think of as a boundary of St. Define nf. to be a unit vector normal to Sf° 



The extrinsic curvature of in St is 

v 

where g% t is the determinant of the metric on St (|A.l| ). The energy of the black hole is 
defined as an integral over the boundary of 

E = lim -— ^ / K bh ^\dS? (A.2) 

r m ax— >oo <57TCt10 J 

dS? = e 2 *R\N c a'fl 2 dSlzdY b \ r=rmax 



Using 



the energy would be 



E = -^RHN c o/f 2 f(r)r ll2 d r (e 2 ^)\ r=rmax J dQ 3 dY 5 

The above expression is divergent and needs to be regularized. Using the background 
renormalization method, we subtract the analogous contribution of a reference background 

Eren = [ " K re f) V\^\dS? (A.3) 

ovrGio J 

For this procedure to be well-defined, the reference background needs to have the same 
asymptotics as the black hole solution, and in particular its metric and fields need to ap- 
proach each other sufficiently fast as r max — * oo. The natural choice for these backgrounds 
is 

ds 2 = e% [-dt 2 + R 2 dQl + -/dr 2 + N c a'dYg] (A.4) 
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where the dilaton and three-form fields are exactly the same as for the black hole solutions 
( |2.1| ), The radius R of the three-sphere is also determined by the 3- form flux ( [2.2[ ), The 
metrics ( |A.4[ ) are also solutions of the equations of motion of the type IIB action coupled 
to the flavor DBI action. 

First of all, the Euclidean time direction of the reference solution has to be compactified 
and its radius at r = r max has to match that of the black hole solution 

1 9tt, ref (j'max ) I — g^bh/ f(r max *) . 

Then, for the metric and fields of the reference background to match those of the black 
hole solution we set 

<f>ref(r) = 4>bh(r) , Rref = Rbh ■ 

As we did above, we can evaluate the extrinsic curvature of a boundary surface Sf° in 
the reference background 

Kref=T l/2 e-^d r {e^)\ r=rmax 
The renormalized energy becomes 



Eren = — „ [ {Kbh ~ K re t) ^\g tt \dS\ 

»7TGio J 



'00 



lim ^r 1/2 R\N c a')^ f{rf' 2 {f{rf/ 2 -l)d r {e 2 % 



dn 3 dY 5 



= -J—e^+^^^R^Nca'f/ 2 [ dQ 3 dY 5 (A.5) 
8vrGi J 

where <f>o + tq is the value of the dilaton at the horizon radius Tq. 

The dependence of E ren on the temperature comes solely from the fact that the radius R 
and 7 are functions of T 

Eren = 8^ e20 ° +2rO(Arca ' )5/27_1/2(T)i?3(r) / dU3dY5 

2<t> +2r ( A2 T 

^ >-N?V(S 3 )V(Y 5 )—, (A.6) 



8vrGio c T H ' 

where we have defined V{S 3 ) = R 3 f dfi 3 and V(Y 5 ) = J dY 5 . 

The energy of the black hole solutions is proportional to their horizon area 

A = e 2<t>0+2ro {N c a'f /2 V{S 3 ) V(Y 5 ) , (A.7) 

and the free energy of the black holes is vanishing 

F = E ren — TS = , 

where the entropy is 

S = ^ 
4G 10 ' 

The fact that the free energy vanishes can also be derived directly from a computation of 
the renormalized Euclidean action of the black holes Z, that is related to the free energy 
by 1 = FjT. 
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B. Equations of motion for metric and dilaton fluctuations 

The equations of motion for the metric and dilaton fluctuations around the black hole 
background can be deduced from the following five-dimensional action 



/ 



* -4/3 



The background black hole metric and dilaton are 

-{l-u)dt 2 + dxl + N c a' 



ds 2 =e^' 



du 2 



4u 2 (l - u) 



4>b = ~ ^ogu 
The equations of motion are 



1 - 4„ ,„ , 2 _ 4/3( 



Rpu - -^g^R - -d^4>d u 4> + -(d(f>) g»„ 



N r 



9txv = 



V 2 0-_ e - 4 /3 =0 

We will consider fluctuations Sg^ = h^^u^t, z), S(j)(u,t,z) in the form 

ha =e -^*+^ e 4/3^( 1 _ u )H tt {u) 



e- iu)t+igz e 4/3<t>b H zz {u) 
60 = e - iu}t+i i z ip(u) 
h tx =h ty = e- wt+l « z e A l^ b H tx {u) 



h tz 
h zz 



h 



xy 



D -iuit+iqz 4/3 



H xy (u) 



and the following choice of gauge 

htu — h xu — /i. 



yu 



h. 







(B.l) 



(B.2) 

(B.3) 
(B.4) 



(B.5) 



(B.6) 



The sound mode metric and dilaton fluctuations {H tt , H tz , H zz , H aa , ip} satisfy the following 
system of linear differential equations 



N c qu) 



4u 2 (l - 



-H„ 







u 



-i2u(i - «)«, + K z ) + Q<Ka + Kz) + 16(1 - «y + ^^H aa -% = q 

u u 
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12u 2 (l - u) 2 (K a - 2H» + 2H" ZZ ) - 4u(l - u){3uH' aa - 9uH' tt + QuH' zz + 8(1 - u)<f/) 
+3N c (uj 2 - q 2 (l - u))H aa + 6iV c g 2 (l - u)H u + 12N c qcvH tz + QN c uj 2 H zz + 32(1 - u)<p = 



12(1 - u) 2 {K a - ~ 12(1 - u)H' aa + 18(1 - u)H' tt - 16(1 " } V 

u 

, 3iV^ 2 16(1 -n) 

H 5 -Haa H 5 ¥> = 

24(l-n) 2 /-24(l- U )^'-^ L( H > H > +H > )+ ( 8 (i _ u ) + _ (1 _ u) 2) ) Q 

6w(l - u)(H' aa + fl£j + 6?(1 - u)H' tz + 3w(F aa + H zz ) + 6g# t2 - 8 ^ (1 - U) if = 

u 

Qq{\ - u)(H' aa - H' tt ) - 6coH' tz + 3qH tt - ^ l ~ U \ = 



2u(l - u) ( 3(u - 2)(H' aa + ff^) + 6(1 - u)H' tt + 8(1 - u)y/ ) 
+3^ c (w 2 -(l-n)g 2 )F aa +3^ 2 (l-u)^i+6^ c gw^+3^ 2 ^ + 16(l-n)<^ = (B.7) 
Conversely, the shear mode fluctuations {Hf X , H zx } satisfy the system 

uH' tx + ^(1 - u)H' zx = , 

^» " 4(1 - «) M > g g » + ^ g *> = ' 
H" zx - ^-H< zx + ^ C _ U? (^H ZX + uqH tx ) = . (B.8) 
Finally, the scalar mode fluctuation H xy satisfies the equation 

The functions {H tt , H tz , H zz , H tx , H zx , H xy , ip] are in general not invariant under the resid- 
ual diffeomorphisms that preserve the gauge choice (|B.6|) 



x» -» x» + e , 
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It is convenient to introduce the following gauge invariant variables [28] 

Zsound = q 2 H a + 2quH tz + cu 2 H zz + q 2 (l -u)(l+ - - f ) ^ , (B.10) 

\ 2(1 — u) q z {l — u)/ 2 

Z 4> = ip--H aa , (B.ll) 

o 

Z shear = lH tx + U>H ZX , (B.12) 

Z scalar = H X y , (B.13) 

The differential systems then lead to 

I , N c (uj 2 - (1 - u)q 2 ) 



q 2 u - 2uj 2 , iV c (4g 4 - 8q 2 uJ 2 ) 

(l-u)(2u 2 + q 2 (u-2)) sound ~ 3(l-u)(2u 2 +q 2 (u-2)) 



4g V (1 -u)- N c {q\u - 2){u - 1) + q 2 u 2 (3u - 4) + 2u 4 ) 

4u 2 (l-u) 2 (2u 2 + q 2 (u-2)) sound ~ ' [ ' 

7 „ ^ 7' jV c (^ 2 -(l-^)g 2 ) ^ _ 

^hear (1 _ n)(w 2_ (? 2 (1 _ n)) Z S W+ 4u 2 (1 _ n) 2 ^ W - U . (tf.lb) 

^scalar ^ _ ^ scalar ' 4-u 2 (l — u) 2 scalar ~ u ■ \ D - LI J 

Note that and 2J sca ; ar satisfy the same differential equation, which corresponds to the 
equation of a minimally coupled scalar in the black hole background. 
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